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Abstract 

For a multicomponent wave field propagating into a multidimensional con- 
version region, the rays are shown to be helical, in general. For a ray-based 
quantity to have a fundamental physical meaning it must be invariant under 
two groups of transformations: congruence transformations (which shuffle 
components of the multi-component wave field) and canonical transforma- 
tions (which act on the ray phase space). It is shown that for conversion 
between two waves there is a new invariant not previously discussed: the 
intrinsic helicity of the ray. 
PACS numbers: 03.65 Sq,42.15 Dp,42.25 Bs 
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For linear propagation of an iV-component wave field in a (weakly) non-uniform multi- 
dimensional medium, the field is usefully represented by a family of rays, along which the 
wave phase, amplitude, and polarization propagate. (The theory is referred to as ray-tracing, 
WKB, or eikonal.) In a local region where two waves of different polarization have (nearly) 
the same frequency and wave- vector, resonant conversion occurs, and a ray of one wave splits 
into two rays, one of each type. In a conversion region, the WKB theory breaks down and a 
local theory must be developed which governs the pairwise interaction of the two resonant 
waves. Such a theory is represented by a 2 x 2 local wave equation. 

The ray equations are Hamiltonian with the determinant of the full NxN dispersion ma- 
trix playing the role of the Hamiltonian. In a conversion region, the local ray geometry can 
be used to guide the reduction from the original NxN theory to the local 2x2 form [1,2], 
construct the local transformations which recasts the local 2x2 wave equation into the 
simplest possible form, and find the local coordinates in which the 2x2 wave equation sepa- 
rates. The new results reported here are: 1) In a conversion region, the six-dimensional ray 
phase space is locally foliated by four-dimensional invariant subspaces. Within these four- 
dimensional subspaces, the ray dynamics is hyperbolic in one two-dimensional subspace, and 
elliptic in the other (i.e. pure hyperbolic motion is not possible). 2) There is one quantity, 
not previously identified, that is invariant under all local congruence transformations and 
canonical transformations: the ray helicity. 

Linear wave conversion is ubiquitous throughout physics (see, for example, references 
given in [2,3]). While there is a large literature on conversion in one dimension, there are 
relatively few studies of the multi-dimensional case (see Refs. [4-14] in physics, and [15-17] 
in mathematics). 

We start with a general iV-component field ^ = (ipi,ip2, ■ ■ ■ , ^n) in a weakly non-uniform 
medium. Assume the medium has three spatial dimensions for concreteness and is time- 
stationary, but note that the methods described can be extended to include non- stationary 
media as well [18]. Denote a point in space as x = (xi, x 2 , x 3 ). The wave equation can 
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always be cast into the standard form [6]: 

J d 3 x'dt'D mn (x, x', t - t')^ n {x.', t') = 0, m — 1,2 ... N. (1) 

We assume that the wave kernel D is an N x N matrix-valued fuction of its arguments 
satisfying D mn (x, x', t—t') = D* m (x', x, t'—t), which gives non-dissipative wave propagation. 
The goal is to find the multi-component field ^(x, t) throughout a given space-time region 
with fitting to appropriate initial/boundary conditions. A standard tool for this analysis is 
the eikonal, or WKB, method. An ansatz is used of the form \I/(x, t) = e l6, ^~^*^(x)e(x), 
where 6*(x) is a rapidly varying phase, V»(x) is a slowly varying scalar amplitude function 
(assumed real), and e(x) is a slowly varying (complex) polarization vector. The WKB 
approximation breaks down near caustics and in conversion regions. The theory of caustics 
is well-developed [19,20], while the conversion problem in multi-dimensions is our present 
concern. 

Starting with the wave kernel and using the Weyl calculus [12,18,21,22], we next con- 
struct the dispersion matrix, an N x N matrix-valued function on the ray phase space 
denoted D(x, k). From now on we denote a point in the six-dimensional ray phase space 
by z = (x, k) = (zi, Z2, Z3, Z4, z$, zq). If (1) is conservative (as we assume), then D(z) is 
Hermitian at each point z: D^(z) = D(z). Ray propagation requires D(z) = det(D) = 0. 
This single scalar condition upon the six coordinates z defines the dispersion surface which 
is, generically, a smooth five- dimensional surface, though it may have local singularities. 
Conversion occurs in the vicinity of the spatial point x* where two different WKB waves 
of frequency cu, with distinct polarization and dispersion characteristics, have nearly equal 
wavevectors [2]. In the ray phase space, this requires that two branches of the dispersion 
surface are in close proximity near the point z* = (x*, k*). 

Away from conversion regions (and caustics), the amplitude and polarization vary slowly 
following a ray. Within conversion regions, however, the polarization and amplitude vary 
rapidly and the WKB ansatz is no longer valid. Instead, the field locally has the form 
^(x, t) = e ik *-( x - x *)-^* [0 Q (x)e Q + ^(x)^], where e a and eg are the uncoupled polariza- 
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tions evaluated at the conversion point. The algorithm for finding these two constant polar- 
ization vectors has been discussed elsewhere (see [6,8] for a proposed algorithm). We assume 
the uncoupled polarizations are given. The complex scalar amplitude functions, ip a and ipp, 
include all the effects of the coupling and have rapid variation in both amplitude and phase 
in the conversion region, but (after multiplication by exp[ik* • (x — x*) — iut]) they connect 
smoothly onto the incoming and outgoing WKB wavef unctions. 

The uncoupled polarizations are used to reduce the full N x N dispersion matrix to the 
following 2x2 reduced dispersion matrix 



D(z) 



^ D aa (z) D a p(z) ^ 



V 



(2) 



/ 



D* af3 (z) D m {z) 

with Djk(z) = e] • D(z) ■ e& j,k = (a, (3). We now consider only the reduced dispersion 
matrix (which is also Hermitian). Following Littlejohn and Flynn [3,11,23], we use the fact 
that any 2x2 Hermitian matrix can be expanded using the Pauli matrices as a basis: 



D(z) 



iyzK 



/ B (z) + B 3 (z) B 1 (z) + iB 2 (z)^ 



V 



(3) 



/ 



B 1 (z)-iB 2 (z) B (z)-B 3 (z 
Here, the components of the 'four- vector' B = (B , Bi, B 2 , B 3 ) are real scalar functions of z 
which are assumed to be independent in the region of interest. Taking the determinant gives 
det(D) = Bl-B\-B\-Bl = rfBpB,,, with the Minkowski tensor 77 = diag(l, -1, -1, -1). 
Note that det(D) = implies that B must lie on the 'light' cone in '5-space'. 

Rays are propagated on the dispersion surface using D(z) = B^B^z) as the ray 
Hamiltonian. Hamilton's equations are most compactly written in terms of the Poisson 
bracket. For any two scalar functions, /(z) and g(z), the Poisson bracket is defined to be 
{/, g} = V z f ■ J ■ V z g where the 6x6 matrix J is defined as 

f A 

-1 

with '0' and '1' the 3x3 null and identity matrices, respectively. The ray Hamiltonian D(z) 
generates the ray evolution equations via: 



J 



V 



(4) 



dz 

da 



{D,z} = -J-V Z D. 



(5) 



Thus, x 



Vfc-D and k = V X D. Following a ray, any scalar function /(z) changes as 



/ = {D, /} = {i] up B u B p , /}. In particular, the coordinates in 'B-space' change as 




/i = 0,l,2,3. 



(6) 



Here 0? = r] up Vt pp is a 4 x 4 matrix composed of all pairwise Poisson brackets of the compo- 
nents -B M (z): Q pfl (z) = {-Bp, B p } = —{B^, B p }. The matrix f2 plays a fundamental role in the 
theory; since all of its entries are Poisson brackets, the entire matrix is invariant under canon- 
ical transformations. Infinitesimal Lorentz transformations are generated by anti-symmetric 



Hence, we can associate a one-parameter family of Lorentz transformations with the ray 
propagation. The key idea is that the equations (6) are simply Hamilton's equations writ- 
ten in a non-canonical coordinate system. (Recalling that phase space is six-dimensional, 
and there are only four B pi we must supplement (6) by two further evolution equations for 
another pair of coordinates, as will be shown.) Note also that, when B — > B' = AB for any 
constant linear transformation A (not just a Lorentz transformation), the matrix of Poisson 
brackets Q transforms via Q' = AVtA. 

The geometrical picture is as follows (see Figure (1)): the ray phase space is six- 
dimensional and plays the role of the base space. Over each point z there is a fiber consisting 
of the space of 2 x 2 Hermitian matrices ('D-space'). Equivalently, the fiber consists of the 
space of four- vectors B ('.B-space'). The fiber space is four-dimensional. For each z, assign 
a particular D(z) and assume that this assignment changes smoothly as we vary the base 
point z. This assignment defines a six-dimensional surface, denoted S. We note that a ray 
trajectory in the phase space (a solution of (5)) is associated with a well-defined curve in S 
and, hence, with a well-defined curve in the fiber space. This curve obeys (6) in .B-space. 

The mapping between the six-dimensional section S and the six- dimensional ray phase 
space is smooth and one-to-one. However, the mapping from S to the four- dimensional D- 



matrices (meaning that the 4x4 matrix A(cr) = 1 + ar/il satisfies A(a)r]A(a) = rj + 0(a 2 )). 
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and .B-spaces cannot be one-to-one. This can be clarified by an appropriate choice of local 
coordinates. We have assumed that the four components B^(z) are locally independent func- 
tions of z, hence we can use them as four local (non-canonical) coordinates in the ray phase 
space. These are supplemented by two further independent coordinates, call them i*\(z) and 
F 2 (z). We can choose F l and F 2 to satisfy {F 1: B^} = {F 2 , B^} = 0, fj, = 0, . . . 3, but Fi 
and F 2 are otherwise arbitrary (note that this implies {F 1 ,F 2 } ^ 0). This gives us a local 
six-dimensional (non-canonical) coordinate system (B (z), -Bi(z), B 2 (z), B 3 (z), Fi(z), F 2 (z)). 
The ^-coordinates change via (6), while the F 1 and F 2 do not change along the ray since 
they have zero Poisson bracket with the B^s. This implies that the rays lie in surfaces of 
constant (F 1 ,F 2 ) and, therefore, there is a natural foliation of the six-dimensional phase 
space, with each four- dimensional leaf labeled by the two invariants F 1 and F 2 . Within each 
four-dimensional leaf, the three-dimensional image of the light cone B^B^(z) = is a slice 
of the dispersion surface D(z) = 0. Note that we do not need to assume that we are near 
the apex of this cone, nor have we linearized the coordinate functions. These properties of 
the local ray dynamics are purely a consequence of the fact that the ray Hamiltonian is the 
determinant of a generic 2x2 Hermitian matrix. 

We now fix attention upon a particular leaf (i.e. fixed values of (Fi, F 2 )). After restricting 
to a fixed leaf, the ray motion can either be viewed in the phase space (where the ray 
equations (5) are canonical), or in £?-space (where the ray equations (6) are non-canonical). 
Extracting the geometrical invariants, and their physical implications, is more direct using 
the non-canonical coordinates. 

Under a congruence transformation Q with constant entries, the 2x2 reduced dispersion 
matrix transforms via D -> D' = Qt D Q, and the determinant as def(D') = Q 2 det(T)) 
where det(Q) = Q. First consider Q with unit determinant, Q — 1. Direct calculation shows 
that the related four-vector B and matrix Q transform to B' u = A^B^ and f2' = Af2A, re- 
spectively, where A is the 4x4 Lorentz matrix with entries [11,21] (A -1 )^ = \tr (o^QV^Q) . 
Notice that two Q's that differ by an overall minus sign are related to the same Lorentz 
transformation. 
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The Lorentz transformations leave the Minkowski tensor r] invariant Ar/A = r\. The 
Minkowski tensor satisfies rf = 1, implying that A~ x = rjAi]. Therefore, 7/ft transforms 
via the similarity transformation i]Q f = A -1 (77ft) A, and its characteristic polynomial is 
invariant under Lorentz transformations. Some algebra shows that P(X) = det(r)Q — A) = 
A 4 — |tr((?7f2) 2 )A 2 + det(rjQ). If we use the standard parametrization for generators of the 
Lorentz group [24]: 



ft 



(7) 



71 72 73 
-71 —u 3 uj 2 

-72 U> 3 -Wi 

^ -73 -uj 2 Wi j 

then we find P(X) = A 4 + (uj 2 — 7 2 )A 2 — (7 • uj) 2 . Notice that P(X) depends only upon A 2 , 
hence the roots of P(X) = come in plus-minus pairs. Generically, P(0) is strictly negative, 
while -P(A 2 ) is concave upward as a function of A 2 . Therefore, as a function of A 2 , P = 
will have one negative root and one positive root. This implies, in turn, that as a function 
of A, P = will have a pair of ± pure imaginary roots, and a pair of ± real roots. This 
implies that the ray motion will be a combination of elliptic and hyperbolic behaviors. 

This parametrization (7) provides a useful physical interpretation for the meaning of ft. 
Fix a point z = zo on the ray where fto = ft(zo) is evaluated. Now use fto to generate a 
one-parameter family of Lorentz transformations in £>-space (a is the ray orbit parameter) 
B(a) = A(a)B(0) = exp(arjil )B(0). This provides a local approximation to the ray orbit 
in S-space. In that space, the ray evolution will be a combination of a rotation (generated 
by the u;-subspace) and a 'boost', generated by the 7-subspace [24]. Hence, the orbit in 
5-space will generically be a combination of elliptic and hyperbolic motions. But, a ray 
orbit in 5-space is the smooth image of a ray orbit in the original phase space, hence ray 
orbits generated by generic 2x2 dispersion matrices will have a combination of elliptic and 
hyperbolic motions. Purely hyperbolic motion is not possible. 

The above discussion demonstrates that ray dynamics in a multi-dimensional conversion 
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is a combination of elliptic and hyperbolic motions (and degenerate versions of these motions 
such as occurs, for example, when uj — 0). However, we have not yet exhausted all possible 
congruence transformations. Now consider congruence transformations that are pure scale 
transformations (i.e. diagonal matrices with |rfet(Q)| =Q^1). From det(D') = Q 2 det(D) 
we have rf v B'^B^ = Q 2 rf v B il B v . Thus the light-cone is invariant, but the numerical value of 
the determinant off the light-cone can change. Using B' = AB we find that A now preserves 
the Minkowski tensor only to an overall scale factor: Ar/A = Q 2 rj, implying A is not a 
Lorentz transformation, but a conformal one. We still have Q' = AQA, and therefore r]Q' = 
Q 2 A~ 1 r]QA. This implies that the characteristic polynomial of i]Q' is P'(X) = det(r]Q! — A) = 
det{A~\Q 2 r]n - \)A) = Q 8 P(\/Q 2 ). Therefore P'(X) = A 4 + (to 2 - 7 2 )A 2 Q 4 - (7 • uj) 2 Q 8 , 
and we arrive, finally, at the result that the ratio 

K^^± (8) 
7 • UJ 

is the sole quantity that is invariant under all (constant) congruence transformations of the 
2x2 reduced dispersion matrix. As mentioned earlier, the entries of Q are Poisson brackets, 
hence the entire matrix is invariant under canonical transformations, and we have uncovered 
a new quantity that is invariant under both sets of transformations. 

In a separate paper we shall discuss how to exploit the connection with the Lorentz 
group to construct the normal form of the reduced dispersion matrix. If 7 and uj are not 
already parallel, it is possible to make them so by performing a 'boost' in the direction uj x 7 
with an appropriate choice of boost parameter. Hence, there is a set of frames where uj 
and 7 are parallel. In such a frame, the invariant we have uncovered reduces to a simple 
function K = k — k~ x with k = cu'/^f'. The primes indicate that these are measured in these 
special frames. This reveals that the invariant k has a natural interpretation as the helicity 
of a ray and measures the rate at which it 'rotates' relative to the rate of exponentiation. 
The sign of k (the 'handedness' of the ray helix, as defined in the four-dimensional B- 
space) is also invariant. We shall also discuss how to solve the related 2x2 wave equation 
using generalizations of the fourier transform, and how to incorporate these new results into 
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numerical ray tracing codes. 
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FIG. 1. A conceptual figure showing the relationship between the six-dimensional ray phase 
(represented by the plane at the bottom of the figure), the four-dimensional fiber spaces (repre- 
sented by the l D' and '-B' planes at the left), and the six-dimensional surface D(z) = S. Note that 
a ray in the phase space defines a smooth curve in S and in the fiber spaces. 
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